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The evolution equation for active and sterile neutrinos propagating in general anisotropic or polarized 
background environment is found and solved for a special case when heavy neutrinos do not decouple, 
resulting in non-unitary mixing among light neutrino states. Then new CP violating neutrino 
oscillation effects appear. In contrast to the standard unitary neutrino oscillations these effects 
can be visible even for two flavour neutrino transitions and even if one of the elements of the 
neutrino mixing matrix is equal to zero. They do not necessarily vanish with Sm 2 — > and they are 
different for various pairs of flavour neutrino transitions [y e — » v^), (y^ — * u T ), (y r — > v e ). Neutrino 
oscillations in vacuum and Earth's matter are calculated for some fixed baseline experiments and 
a comparison between unitary and non-unitary oscillations are presented. It is shown, taking into 
account the present experimental constraints, that heavy neutrino states can affect CP and T 
asymmetries. This is especially true in the case of — > v T oscillations. 



I. INTRODUCTION 



In the last few years a common effort of experimentalists and theoreticians yield new informations on the neutrino 
sector. Nowadays we can say that masses and mixing of the three light neutrinos are quite well established. The 
unknown issue is the total number of neutrino species Q. Moreover, additional neutrinos can be both light (with 
masses of the order of electronovolts) or heavy (with masses greater than the Z boson mass) 0. In both cases 
their couplings to the ordinary matter must be much smaller than the couplings of the three known neutrinos. 
Additional light neutrinos are called sterile. Such neutrinos are still permissible from the point of view of primordial 
nucleosynthesis || and seem to be necessary if LSND results || are confirmed by MiniBOONE [[!]. As far as heavy 
neutrinos are concerned, their number is not established. For cosmological reasons they must be unstable || or they 
do not exist at all besides the heavy compact objects with the most exotic physics fij. 

Let U v be the full neutrino mixing matrix and let the submatrix IA (of dimensions (3 + n s ) x (3 + n s )) constitute the 
mixing matrix of light known neutrino states and let the submatrix V (of dimensions (3 + n s ) x iir) be responsible 
for the mixing of the light neutrinos with tir additional heavy states. Then 



Many new physics models predict heavy neutrinos. In case they are very heavy objects at the unification scale, the 
see-saw mechanism explains why the known neutrinos are light and the matrix U, which enters neutrino oscillation 
formula, is practically unitary (light-heavy neutrino mixing representing by the V submatrix is negligible). Without 
the see-saw mechanism the submatrix V could have substantial elements. This situation might be realized with TeV 
neutrino masses. However, even then the construction of the matrix U v is not trivial [g, unless all elements of the 
neutrino mass matrix are of the same order of magnitude and symmetries disconnect light-heavy mixing from the 
ratio of their masses [[| . 

Nowadays, experimental data constraints elements of V by terms of the following form |j| 

c aP = (vvt) Q/3 = J2 (2) 

i— heavy 

so the submatrix U, which decides about the light neutrino oscillation effects is not unitary, and nonunitarity effects 
are conjugate with c a p 

{Urf) ap =5 aP -c aP . (3) 

This means that various channels of neutrino oscillations will depend on the c a p nonunitarity parameters. Some 
of such studies have already been done || |To) . Here we would like to concentrate on the influence of the non unitary 
matrix U on CP-violating neutrino oscillations. In the previous paper [ pd] | we have found that CP violating oscillations 
in vacuum can be affected by c a p parameters. We know that, for conventional unitary neutrino oscillations, the CP 
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violation can occur for three (or larger) number of mixed neutrinos, and is small (vanishes) if any single element of 
the U matrix is small (vanishes). The CP violating effects are govern by one (three) Jarlskog invar iant(s) @ for 
three (four) neutrinos. Finally, they vanish for short baseline oscillations E3| . All such limitations can be avoided in 
the case of nonunitary neutrino oscillations. 

In reality, the influence of matter effects on neutrino oscillations is important. It can happen that the matter effects 
mimic or even screen CP effects [0 . We will also investigate the role of the non- unitary U matrix in the neutrino 
oscillations in matter. 

The paper is organized as follows. First, in Chapter II the equation of motion for neutrino states in medium is 
derived. Such an equation is well known in the case of the unitary neutrino mixing. To the best of our knowledge 
the same equation has never really been found for the non-unitary mixing. Next, we use this equation to find various 
CP and T asymmetries, at first analytically (Chapter III) and then numerically (Chapter IV). Conclusions close the 
paper. 



II. PROPAGATION OF STATES IN MATTER 



It is well known that the interference between scattered and unscattered neutrinos can be crucial for their prop- 
agation in matter, even if the probability of incoherent neutrino scattering is negligibly small. There are several 



derivations of the evolution equation for neutrinos in matter | 
which describes the averaged, coherent neutrino interactions wit 



15fl . In all cases, at first an effective potential V e ff, 
i all background particles, has to be calculated. Using 



this potential, the Dirac's equation for a neutrino bispinor wave function \P can be written 

(i^fy - ro - V o// ) ¥ = 0. (4) 

Taking into account that during the evolution of the flavour states, particle-antiparticle mixing is negligible, neutri- 
nos do not change their spin projection and that they are relativistic particles, a simpler, Schrodinger like evolution 
equation can be found 



ij t \u(t))=H eff \u(t)). (5) 

Usually the effective potential V e ff is calculated in the neutrino flavour basis. This is the traditional approach to 
the three active neutrino mixing. If sterile and/or heavy neutrinos exist, it appears to be more natural to calculate 
V e ff in the eigenmass basis. There are several reasons why it is so. At first, it is not clear conceptually how to 
define, in a consistent way, creation and annihilation operators for flavour states [ fL6| . At second, as a matter of fact, 
there is no such an object as a flavour eigenstate. For quarks, for instance, only eigenmass states are used. In this 
context the only difference between quarks and neutrinos are much smaller 8m 2 's. However, as we will see later, 
the most important thing is that: using the eigenmass basis we will be able to avoid the non-hermitian evolution 
of neutrino flavour states affected by the heavy neutrino sector in a matter. The equation of motion in Chapter III 
can be described by a hermitian Hamiltonian from which real effective neutrino masses follow. Final probabilities of 
flavour changing are affected by the heavy neutrino sector through initial conditions and then nonunitary effective 
neutrino mixing appears. 

In order to find the effective potential V e f / , let us assume neutrino interactions in a general form 



Lcc = 



2\/2sinevi/ " 



J2 E ( X - 76) VvUniW- + h.c, 



(0) 



and 



Lnc = 4 ^ ft e m ,fl \ £ ^ 0- ~ 7s) SUsniZp + 2 J2 [Tzt (1 - 7s) - 2Q f sin 2 S w ] V f Z» \ , (7) 

4 Sin <Oy/ COS <0\y r-r^-, r — 

^*)J=1 f=e,p,n ) 

where n is the number of light (3 + n s ) and heavy (hr) neutrinos (n = 3 + n s + riR), and = J2 a=e „ T {Uu)* ai (U v ) a y 
The coherent neutrino scattering is described in general by three types of Feynman diagrams presented in Fig. [h The 
Higgs' particles exchange diagrams do not contribute as neutrinos are relativistic particles. In the normal matter all 
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diagrams contribute to the neutrino-electron scattering m + e~ — > n k + e~ (/ = e), yet only diagram (a) contributes 
to neutrino-nucleon scattering rij + / — ► + f (/ = p, n). 

At low energies (q 2 <C M|y, M§), the effective interaction of light neutrinos with a background particle f (f — e,p, n) 
(Fig. ^]), after the appropriate Fierz rearrangement, can be written in the form 



r i 3+n 



H( nt {x) = -§ Y. E ["* r -"*] [*/ r ° ($ + 3/a75) */] , (8) 

^ i,fc=l a=V. A 

where Ty(A) = 7^(7^7s)- The couplings and can be calculated from Eqs. |(|f?] and for electrons (/ = e) and 
nucleons (/ = p, n) they are given by 

: U* k Uei + P^ki y~ + 2 sin 2 <d w 
- ptt k i (T 3 f - 2Q f sin 2 6^) , 

-pQkiT 3f , (9) 



9eV — 


~9eA 


9eV — 


ki 
9 e A 


n ki _ 
9fV ~ 


9fA 


9fV ~ 


ki 

~9eA 



where 



, T 3p = -T 3n = 1/2, Q P = 1, Q n = 0. (10) 



p M§ cos 2 @ w 

The global effect of matter - light neutrino interaction can be described by the Hamiltonian 



3+n s 



H int (x) = Y n kV k irii, (11) 



i,k=l 

where 



and 

(K/) M = ^|E/^/(P,A)(M/) H . (13) 

(Af/) fei is the part of the matrix element of the scattering amplitude n% + / — > + / connected with the fermion / 
in the case when particles' momenta and spins are untouched 

(Mi) k . = (f,p, A|# f T a (g% + gf alb ) * f \f,p, A). (14) 

In Eq. [l^, Pf(p, A) is the distribution function for the background particles of spin A and momentum p, normalized 
in such a way that Nf, defined as 

A v ' 

is the number of fermions / in a unit volume (V — 1). Hence, the amplitude (M/j, . must be calculated for a single 
fermion in V — 1 (then for bispinors u we have u'u — 1) and |T?j] 
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( M C - - ( M C = & (I) + m ^ (I) ' < i6 > 

where Ef,mf and Sj — (^p\,\mf + J denote the fermion's / energy, mass and spin four- vector, respec- 

tively. The obtained relation between the vector My and the axial-vector Ma amplitudes is the consequence of the 
V-A form of interactions (Eqs. |^J^)- Now, using Eqs. |l3| and [l6] we can write the effective potential V ki (Eq. |l2|) 



hi 



/. / / ^ f_ \ | ry-r> „ / ^ f ' 



9fv{^-)+m f g fv , E 



(17) 



where the average value < z > defined by 



< z >= ^ E / ^3^® X) (18) 

describes quantities averaged over the fermion distributions Pf(p, A). In this way the set of coupled Dirac's equations 
for all light neutrinos propagating in the matter is obtained 

3+n s 

2J (i^dfihi ~ m k S ki - V k i) Ui = 0, fc = 1, 3 + n s . (19) 

i=l 

As we describe the propagation of light relativistic neutrinos only, a simpler Schrodinger-like evolution equation 
can be found. Assuming that k 2 ~ _E 2 >> mf ~ mi\Vki \ — \Vki\ 2 , we can get, in the momentum representation, the 
Schrodinger like equation for each left-handed (or each right-handed) components of neutrino bispinors 

, 3+n B 

i-* k (M) = H%f(t)9i(k,t). (20) 

i— 1 

The wave function ^-(A;, t) is the neutrino (antincutrino) state \^(t) > with momentum k and hclicity A = — 1(+1), 
written in the eigenmass basis \vj >, =< Vj\^(t) >. The effective Hamiltonian (we will assume from now on that 



all neutrinos have the same momentum but different energies Ej = y k 2 + m 2 ) is equal to 



where 



H e J/ =(k + ^)6 kl +Hir , (21) 

A^Uk^f ^PbUi for Dirac neutrinos, 

H ki t=<v k\J d 3 x H int (x) \vi > = < -(A^ i )*u k J IJ ,PRUi for Dirac antincutrinos, (22) 

A^UkJuPLUi - (A^yUk^PRUi for Majorana neutrinos, 



- d ^=£/K> fei - 

In the relativistic limit, with the additional assumption k = ki = kf , we have 



5 



u k l t j,PLUi\ x= _ l = ttfc7 ft P R iti| A=+1 = 



= 1,- 



1*1. 



and 



Ufc7 A1 PLW i | A=+1 = Mfc7At-P-R u *|A=-l = (°'° 



(23) 



(24) 



so, finally, we have arrived (A^ = (A^Aki)) to the formula 



Auj — -rrAhi for Dirac and Majorana neutrinos with A = — 1, 

HIT ={ n 11 r. - (25) 
"(■^■fct)* |f[(^ fci -'* ^ or Dirac antineutrinos and Majorana neutrinos with A = +1. 

We can clearly see that the effective interaction Hamiltonian is the same for Dirac and Majorana neutrinos and 
that 



rrint 
particle 



— r tj 



antiparticlel 



(26) 



We would like to stress, however, that these properties are not the general rules. They are satisfied because of the 
V-A type of neutrino interactions in Eqs. ^|,^, as it is in the case of relativistic neutrinos for which the scalar and 
pseudo-scalar terms can be neglected. 

Eq. |2^ gives the most general Hamiltonian for an arbitrary number of (3 + n s ) light relativistic neutrinos propagating 
in any background medium and interacting in the V — A way. In what follows we will concentrate on the case of 
unpolarized (< A/ >= 0), isotropic (< pf >= 0) and electrically neutral (N e — N p ^ N n ) medium. Then we arrive 
to the following Hamiltonian (later we will usually put p — 1) 



N e U: k U ei - \ P N n tt k . 



(27) 



The above Hamiltonian always has (3 + n s ) x (3 + n s ) dimensions, independently if heavy neutrinos exist or not. 
Let us first consider two conventional cases with none and with a single sterile neutrino, when no heavy neutrino 
exists (jir = 0, n = 3 + n s ). Then (in both cases) 



3+n s 



(28) 



If no sterile neutrinos are present [n s = 0), then, after removing the redundant diagonal terms in Eqs. ^l],^], the 
well known Hamiltonian is obtained (a, (3 = e, p, t) 



HaP = ^^E^^) + ^F^eS/Se- (29) 

If a single sterile neutrino is present (n s = 1), then another well known Hamiltonian is obtained (a, [3 = e, /z, r, s) 



H, 



a/3 



u- 



2E 



-U 



V2Gi 



a/3 



NeSaeSffe + ~N n 5 as 5f3 6 



(30) 



From now on, we will only consider cases when there is at least one non-decoupling heavy neutrino present (tir > 1). 

If we want to use the full (3 + n s + tir) eigenmass basis, wc need to expand the (3 + n s ) x (3 + n s ) Hamiltonian 
Hki 3 F > gi ven by Eqs. pl],p7|, to proper (3 + n s + ur) x (3 + n s + ur) dimensions, adding zeros 
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ireff f H^^as given by Eqs. 21,^, if both k, i < (3 + n s ), 
kl \ 0, if any of fc,i > (3 + n s ). 



(31) 



These zeros mean that, for the physics of light neutrinos which we are interested in, the energy and momentum 
conservation do not allow heavy neutrinos to be produced nor detected. If we now define the flavour basis as 



3+n s 



3+n s +n R 
i=3+n s + l 



(32) 



we can express the Hamiltonian Eq. |3l] in the flavour representation. 

If no sterile neutrinos are present (n s = 0), but some heavy neutrinos exist, n — 3 + n^, Ur > 1, then (a, j3 — e, /x, r) 



H a/ 3 = [id ^ J + V2G F N e (5 ae 5p e - c a p(5 ae + 5p e ) + c ae c e p) + -N n (2c a p - £ c ai c 7 p) 



(33) 



If a single sterile neutrino (n s = 1) and some heavy neutrinos exist, n — 4 + tir, ur > 1, then (a, f3 = e, fi, t, s) 



a/3 



(34) 



a/3 



\/2G; 



{^ae<5/3e ~ C a/3 ((5 Qe + ^ e ) + C eQ C e/3 } + -iV n < (Sasfyjs + C 



a| g(2 — <5 as — <5g s ) — £ C Q7 C 7 > . 



In both above cases (Eqs. |3^ , p4| ), we present the Hamiltonian in the light neutrino subspace only. The full Hamil- 
tonian is more complicated and contains parts related to the light-heavy neutrino mixing. In the full basis, it is 
represented by a (3 + n s + ur) x (3 + n s + Ur) dimensional matrix 



H a p 



(35) 



aP 



where H e ^ is the (3 + n s ) x (3 + n s ) dimensional Hamiltonian given by Eqs. |2T , and Eq. |3^ or Eq. |34| is the 
most upper left part of the matrix Eq. |35|. Thus, even though we consider the problem of light neutrinos propagation 
only, in the flavour basis Eq. |3^ we have to deal with (3 + n s + ur) x (3 + n s + tir) dimensional matrices. 

It should be stressed here that, up to now, all Hamiltonians in Eqs. ^lf]35| are represented by hermitian matrices. 
The transformations in Eqs. ^,^2| are unitary as we sum over all n available neutrino eigenmass states. 

If we now take into consideration the fact that, due to kinematical reasons, no heavy mass eigenstates in Eq. |32| 
can experimentally be produced, then the properly normalized states \v a >, which correspond to neutrinos produced 
in real experiments, are 



3+n s 



3+n s 



>= a- 1 J2 KiWi >= E >. 



(36) 



i=l 



i=l 



where X n 



J2i=i s \Mai\ 2 = y/l — c aa and U a i = X^Uai- Such states are not orthogonal 

< v a \i>p >^ 0, a ^(3. 



(37) 



Let us notice |T(], ^T) that in this case the usual notion of flavour neutrinos loose its meaning. For example, a 
neutrino which is created with an electron, and is described by the state \u e >, can produce besides an electron also a 
muon or a tau. It is better then to see active neutrinos as particles which are produced together with charged leptons 
of particular flavours (in some charged current weak decays), rather than particles having their own flavours. 

When we write Eq. [2(1 in the basis of experimentally accessible states \v a >, we get 
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(38) 



The Hamiltonian H a p is given by (X a [3 = X a S a /3) 



H = UHU 1 = — U 



( 

8m 2 21 

6m%i 





\ 



J 



U 1 +V2G F ^X 2 



({N e -^f) 

u 2 







V 





2 u 





(39) 



This matrix, in contrast to the previous hermitian representations (in Eqs. |2l|-p5D, is not hermitian. In practice we 
can solve the basic Eq. cither in the experimental \v a > or in the eigenmass \vi > basis (in both cases we deal with 
Hamiltonians having dimensions (3 + n s ) x (3 + n s )). 

Let us assume that at time t = the state |^(0) >= |£ Q (0) > is produced. In order to find the state \v a {t) >, we 
need to solve Eq. with an initial condition 



< vp\i>a(0) >= (X a Xp) 1 (S a (3-C0 a ). 



In the eigenmass basis we have 



d 3+™s 

i-j^ < Vk\Va(t) >= ^ Hki < v i\Va(t) > 
i—1 



where 



H 



1 

2E 



( 

Sm^ 






Sm 2 



:S1 



'41 



v.- 



o 








N n 
2 






V 



and the initial condition is 



<Vk\v a (0) >=X a l U* k 



2 u 





u. 



(40) 



(41) 



(42) 



(43) 



For neutrino propagation in a uniform density medium we can solve the evolution equation analytically. Then it is 
simpler to use the eigenmass basis in which the effective Hamiltonian H is hermitian. We will follow this approach in 
the next chapter. In the case of a medium with varying density, Eq. ^ or Eq. ^ must usually be solved numerically. 
Then any of them, hermitian Eq. ^l] or non-hermitian Eq. [38], with appropriate initial conditions, from respectively 
Eq. H or Eq. |o|, can be used. 

Finally, we should also remember that, the nonorthogonality of the \v a > states has some impact on theoretically 
calculated cross sections for neutrino production and detection. Let us consider, for example, a charged lepton 
production process v a X — > l~pY . The production amplitude is then 



3+n s 



3+n s 



A{v a X -> l-Y) = A- 1 KiM"iX - IpY) ~ A- 1 £ W m U 0l A(v m=o X - l~ p Y) 



(44) 



The A(v m= oX — > 1^ Y) is simply the Standard Model's amplitude describing the process v@X — > la Y . Thus, the 
production cross section is 



(45) 



a{y a X -> IpY) = X~ 2 \S a (3 - c a0 \ 2 a bM (v p X -> l~Y) 

where a SM (vpX — > IgY) is the Standard Model's cross section for the process v$X — > l~ZY ' . For other processes, 
which we do not describe here (like, for example, neutrino elastic scattering), even more complicated "scaling" factors 
appear. 
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III. OSCILLATIONS OF LIGHT NEUTRINOS WITH CP AND T VIOLATING EFFECTS 



In what follows oscillations of three light neutrinos with n# > 1 non-decoupling heavy neutrinos will be considered. 
If heavy neutrinos decouple then the mixing between flavour v a = (V e , i/„, v T ) and mass = (yi, i>2, ^3) neutrinos is 
described by the 3x3 unitary matrix U, namely, 

3 

\v a >=J2KiWi> ■ (46) 

i=l 

For the matrix U the standard paramctrization is taken 



C12C13 s 12 ci 3 s 13 e~ \ 

-S12C23 - Ci 2 S23Sl3e ! ' 5 C12C23 - si2S 2 3Si3e l<5 S23C13 ■ (47) 
S12S23 - C12C23S136" 5 -Cl2S 23 - Si2C2 3 Sl 3 e lS C23C13 / 

Now, in order to implement effects of heavy neutrinos, as discussed in the Introduction (Eqs. the submatrix V 
must be introduced. As the number is unknown, V will be parametrized in a simplified way, using a single effective 
heavy neutrino state. In this way, the effective light-heavy neutrino mixing can be characterized ]TT| ] by three real 
parameters e e ,e M ,e r , and two new (real) phases xiiX2 (responsible for additional CP and T asymmetry effects) 



/ U el U e2 U e3 \ 

U=[ Upl UyH = 

\ Url U T 2 U T3 J 



V = I e' l ^e„ I . i4M 



Then A Q = \J\ — e 2 a w 1 — e^/2 and the 3x3 matrix ti(e a ), to the order of 0{t 2 a ), has the form 



U = 



U e 2 A e , 



3 e tet^i 



(49) 



UriK -U el e-**" 

J lP i(xi-X: . r , 



V -U.xe^-^haer 



U t2 X t - U e2 e lXl e e e T 

— II r ,p i (xi-X2), f 



U t3 \ t - U e3 e lXl e e e T 



where U a i are the unitary mixing matrix elements as in Eq. 47 
Experimental data restrict the light-heavy neutrino mixing elements er a 



e\ < 0.0054, e 2 ^ < 0.0096, e 2 T < 0.016 , (50) 

altogether with their products 

e e e M < 0.0001, e M e T < 0.01. (51) 

There are no constraints on additional CP breaking phases xi an d X2, so we will assume, as in the case of the 
standard CP phase 6 that any values in the range (0, 2tt) are possible. The matrix Qf-i in Eq. [7] is given by (R 
stands for the single effective heavy neutrino) 



JIm = *K - V&Vja, M = 1,2,3, (52) 

and 

V' m = -U et e e - L/^e** 1 ^ - U Ti e?**e T . (53) 

Now, assuming a constant matter density, the evolution equation Eq. ^l] can be solved analytically. First, the 
hermitian effective Hamiltonian 



H e J f = 7rz( m 1 5 ^ + 2V2G F E ( N e U* ek U e% - ~iV„ 



2E V ' _ " * ~ \- -^-' l ek L/l ei ^^rS^ki ) ) (54) 
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can be diagonalized by a unitary transformation 

H eff = ^WUiag{m})W, (55) 

where m? are (real) eigenvalues of 2EH e ^ and W is a matrix build of eigenfunctions of 2EH e ^ . Then Eq. [ll] takes 
the form 

= ^WUiagim^W^it), (56) 

where (the \v a > states are given by Eq. |36|) 

/ < vx\u a {t) > \ 

* a (t) = <V2\Va{t) > ■ (57) 

V < v 3 \i> a (t) > J 

This equation together with the initial condition Eq. ^ can easily be solved giving 

~ 2 

n(t) = ^(^Jne-^'fTOt),,, , (58) 

i 

and the amplitude A a —,f}(L) for v a — > vp neutrino oscillations in matter, after traveling a distance L, is given by 

3 i%2 

A a ^ p (L) =< 9(3(0)\v a (L = t)>=J2 W 0l W* ai e~ l ^ L . (59) 
The non-unitary neutrino mixing matrix W is defined as 

W ai = (UW*) ai = X-^UakWik = Kx^Wai. (60) 

k 

The final transition probability P a ^p{L) = \A a ^f3(L)\ 2 is the following 

P a ->p(L) = -^-2 I (Safi ~ (W t ) Q ) 2 - 4 ^ i?^ sin 2 A lk + 87^ sin A 21 sin A 3 i sin A 32 

+ 2 [4 1 jsin2A 3 i+4 2 jsin2A 32 ]}, (61) 

where 

= Re [W ai W lSk W: k W; M ] , (62) 
7^ = 7m [W a ^ fc W* fe W£] , (63) 
A$ (e Q ) = 7m [W^c^] , (64) 

and 

A.- 1 267 ( ^-^ )[e ^ ]LM (65) 
A lfc - 1.267 . (65) 

It is interesting to notice that |HJ ], while in the case of unitary mixing matrix (njj = 0) we always have 
E/3=e n t( s) Pol— — lj it is no longer true when tir > 1. In this case, as a consequence of nonorthogonality of 
the \v a > states, this sum can be bigger or smaller than 1, and its value changes with neutrino energy and distance 
(time). 

In agreement with Eq. ^5], the transition probability for Dirac antineutrino or Majorana neutrino with A = +1 can 
be obtained from Eq. |6l| after the replacement 

P a _>g(£) = P a ^{L-U -^W,G F ^ —Gp). (66) 
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Then conventional CP and T violation probability differences are 



AP 



,cp _ p p _ 



a P k i>fe 



R%(G F ) sin 2 A lk (G F ) - R%(-G F ) sin 2 A tk (~G F ) 



Q(G F ) sin A 21 (G F ) sin A 31 (G F ) sin A 32 (G F ) 



+ P a % [—Gp] sin A 21 (-G F ) sin A 31 (-G F ) sin A 32 (— G F ) 

A% (G F ) sin 2A 31 (G F ) + A™ (G F ) sin 2A 32 (G F ) 
- (~Gf) sin 2 A 3 i (-G F ) + A<g (-G F ) sin 2A 32 (-G F ) 



and 



APj_^ = P a ^ (i - P^ a = 1 4712 gin Aai sin Agi gin + A W gin 2Aai + A (2) sin 2 As2 J 



(67) 



(68) 



Nonunitarity of the W matrix produces two types of effects. At first, all conventional quantities such as R 1 ^,!^ 

depend on e a and new CP phases \i- At second, new terms proportional to appear. 

The first effect can mainly be seen in numerical analysis. The presence of the additional term is more spectacular 
and can be analyzed directly. For neutrino oscillations in vacuum (G F = 0) new terms do not change the relation 
between AP^p and APj_^, and they are equal 



AP£J^p (vacuum) = A P^p (vacuum). 



For a = (3, A^l = and 



AP^ a (matter) — AP^ a (vacuum) = AP a _^ a (vacuum) = 



(69) 



(70) 



In the normal medium, which is matter and not antimatter, AP^ a (matter) 7^ 0. Heavy neutrinos will make this 
effect stronger. For long baseline (LBL) neutrino oscillations 



A L bl ~ A31 ~ A 32 ~ O(l), A 2 i~0, 
and, in contrary to the unitary oscillations, 

AP£^p(vacuum) — AP£Ilp(vacuum) 7^ 0, 



APl^Jmatter) 



(A^ d +A^>)sm2A LBL 



t(i) 



1 a/3 



- A TTW Im i W ^s W P3 c* a[j ) sin 2A LBL . 

A a A p 



(71) 
(72) 

(73) 



This means that CP and T asymmetries appear even for two flavour neutrino transitions. Furthermore, for unitary 
3 flavour neutrino oscillations, moduli of all Jarlskog invariants are equal and, as a consequence, all T asymmetries 
(equivalent to CP asymmetries in vacuum) are equal as well 



AP 1 = AP 1 = AP 1 . 



(74) 



In addition, if any element of the mixing matrix is small (vanishes) then the above asymmetries are also small 
(vanish). For V^0 (and as a consequence V 7^ 0), there is 



1% = 1% + Im[W^W; k (WV' J V*W% k ] = If e - Im[W ei W* ek (WV 2 V*W^) ik ], 

lap = i 2 a p-i™iw: 2 w 02 c* af3 ] 



(75) 
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The above relations and terms proportional to A^l in Eq. ^38| imply 



AP e % ^ APj_ ^ AP T ^ e . (76) 



From Eq. [75] follows that even when some of the U matrix elements vanish, the asymmetry AP^^„ can be nonzero. 



For instance, if W e 3 — then I l 3 3 = (i = 1,2, /3 = /z, r), but remaining five CP invariants, where is absent, do 



not vanish. 



IV. NUMERICAL RESULTS 



Here we will present some numerical analysis of the standard and nonstandard CP and T violating effects for real 
future neutrino oscillation experiments. In order to check the effect of non-decoupling heavy neutrinos in neutrino 
oscillations, the CP and T asymmetries for two baselines L=295 km and L—732 km are calculated. Neutrino energy is 
allowed to vary, according to the experimental conditions, in the range E = 0.1-^30 GeV for v e — > and E = 1.78-^30 
GeV for — > v T . These baselines are planed for several future experiments (e.g. JHF and SJHF in Japan with E ~ 1 
GeV 0; ICARUS || and OPERA j2l| in Europe, E ~ 20 GeV; MINOS g in USA, E ~ 10 GeV; and SuperNuMi 
|p3| in USA with P ~ 3,7, 15 GeV). In all these experiments neutrino beams from the pion decays will be used, so 
they are mostly muon neutrino and antineutrino beams. Neutrino factories will give in addition electron neutrino and 



antineutrino beams. In Figs. 2a-6c the probability difference AP^_^ for 
also standard quantities defined as 



and 



is shown. For v u 



\CP 



A 



are presented and matter effects are included 



ap cp 



AP 



P{vp, —> Vr) + P(v T — > Vu 



(77) 
(78) 



A e [eV 2 ] = 2V2G F N e E = 7.63 x 10~ 5 
A n [eV 2 ] = V2G F N n E = 7.63 x 10~ 5 



P 



g/cm 3 
P 



g/cm 3 



Y, 



0.5 

1 1 -Y t 



E 



GeV 
E 
GeV 



(79) 
(80) 



For L=250 km and L=732 km neutrinos pass only the first shell of the Earth's interior pif with a constant density 
p = 2.6 g/cm 3 and Y e = 0.494. Then A e [eV 2 ] = 1.96 x 10" 4 [E/GeV] and A n [eV 2 ] = 1.0 x 10~ 4 [E/GeV]. 
The oscillation parameters are taken from the best LMA fit values |^5j tan 2 O23 = 1.4, Sm 2 2 = 3.1 • 10~ 3 eV 2 , 

3.3 • 10^ 5 eU 2 and tan 2 26i3 = 0.005. 

and e T (Eq. Bl) are small and satisfy the experimental constraints given by Eqs. |5^,pl|. In 



tan 2 6i2 = 0.36, 8m 2 2l = 
The parameters e e , e p 



agreement with the above constraints two sets of e a parameters will be discussed 



{A) : e e ~ 0.001, e n = e T = 0.1, (81) 
(B) : e e =e M ~0.01, e T = 0.1. (82) 

Figs. 2a,2b present AP^^ (L=250 km) for neutrino oscillations in vacuum (Fig. 2a) and matter (Fig. 2b), respec- 
tively. We can see that matter effects are very weak. The hatched regions in these figures, and in all following figures, 
describe the classical unitary neutrino oscillations, with < S < 2ir, and they are similar in both cases. In Fig. 2b this 
region is only slightly asymmetric. As parameters e a which describe non-unitary oscillations are very small, deviation 
of AP^ from unitary oscillations is very weak. The shaded regions in Figs. 2a,2b, and in all following figures, give 
the allowed range of AP^J/ a for the non- unitary case with additional CP breaking phases xi > Xi (Eq. |i^) that change 
in the domain < Xi < 2-7T. It is interesting to see the same CP violating quantity AP^J/^ for 813 — > 0. In vacuum this 
quantity equals to zero. In Fig. 2c a possible range of AP,^^ is depicted for non-unitary oscillations. In agreement 
with our previous discussion (Eq. |7^) such a quantity does not vanish. However, it is small, which comes out of strong 
experimental bounds Eqs. pT],p2l Situation does not change qualitatively in the matter case (Fig. 2d), except that this 
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time unitary oscillations can be nonzero. Such miserable effects cannot be detected in future neutrino experiments. 
Of course, A^ and -A;T effects alone can be large (see e.g. We would like to stress, however, that A^^g 

and A^^p must be discussed together with AP^p quantities. Only then we can say if these effects can really be 
measured in experiments. Such a discussion for — > v T transitions will follow in the next figures. 

In Figs. 3a,3b the APf L Z T and Af t ^ T quantities in vacuum are presented as functions of neutrino energy for L=250 
km and the first set of e a parameters from Eq. |8l]. Here the CP violating effects are quite large. For example, for E=2 
GeV, AP^ T e (-0.0003, +0.0003) in the unitary case (hatched region in Fig. 3a) and AP^ T e (-0.017, +0.017) 
in the non-unitary case (shaded region in Fig. 3a). We do not present results for e Aj^_ tT here as in vacuum A£_^„ = 
A^^g- We can see that for higher neutrino energies non- unitary effects can be very large (and increase with neutrino 
energy). Unfortunately, oscillation probabilities themselves are getting smaller. 

In Figs. 4a, 4b the same quantities as in Figs. 3a,3b are presented but for the second set of e Q parameters (Eq. |32"| ). 
As e M is 10 times smaller, also nonunitary effects are smaller by the same factor. 

For L=250 km the effects of neutrino interactions in the Earth's matter are small and all presented quantities are 
almost the same as in vacuum. Thus we do not present them here. However, for longer baseline experiments (L=732 
km) matter effects can already be seen. In Figs. 5a, 5b and Figs. 6a-c the AP9^ T , A^ T and A^^ T quantities are 
presented for — ► v T transitions in vacuum and matter, respectively. The e a parameters are taken according to 
Eq. ^l| Again, the non-unitary effects are large. For the unitary case, the symmetric region of AP^ T in vacuum 
(Fig. 5a) becomes asymmetric (negative) in matter (Fig. 6a). The full range of Af~^ T for non- unitary oscillations 
shifts also slightly toward negative values (compare the shaded regions in Fig. 5a and Fig. 6a). Once again Aj^ T 
(Fig. 5b and Fig. 6b) and Aj^_ tT (Fig. 6c) asymmetries are very similar to each other. They are larger for higher 
neutrino energies but, unfortunately, not because AP^ is larger, but because the probabilities for (anti)neutrino 
oscillations are getting smaller. 

V. CONCLUSIONS 

The evolution equation for neutrinos propagating in matter is found in the case when heavy neutrinos do not 
decouple and, as a result, the mixing matrix between light neutrinos is not unitary. 

The neutrino propagation is described by a hermitian Hamiltonian found in the eigenmass basis of light neutrinos. 
This basis is the most convenient to calculate physical effects. Two other basis which we consider, that means the full 
orthogonal basis which contains both light and heavy neutrino states, and the "experimental" basis with nonorthogonal 
light neutrinos only, are more complicated when applied to numerics. In the first case, the Hamiltonian is hermitian 
but has full (3 + n s + Ur) x (3 + n s +n^) dimensions. In the second case, it has smaller (3 + n s ) x (3 + n s ) dimensions, 
but is represented by a non-hermitian matrix. 

If heavy neutrinos do not decouple, the notion of neutrino flavour loose its meaning. In such case it is better 
to see active light neutrinos as particles produced together with charged leptons of particular flavours rather then 
particles having their own flavours. Moreover, the nonorthogonality of such neutrino states has also some impact on 
theoretically calculated cross sections for neutrino production and detection. 

We have found the equation which describes the propagation of light neutrinos in any polarized, nonhomogeneous, 
charged matter. Numerical analysis is done in the simpler case of unpolarized, isotropic, and electrically neutral 
matter. 

The non-unitary neutrino mixing is especially important for CP and T violating effects in neutrino oscillation. Two 
additional CP phases, which appear in the light-heavy neutrino mixing matrix, have crucial consequences in such 
phenomena. 

In comparison with normal unitary neutrino oscillation, the CP and T violation are generally larger and appear 
even where standard effects are very small or vanish. The new CP and T violation effects depend on the strength 
of light-heavy neutrino mixing. For present experimental bounds on the non-decoupling parameters, the — ► v r 
oscillation is especially sensitive to the new CP effects. For v e — > oscillation, where the bound on the c efl parameters 
coming from the non observation of the fx —* ej decay are very stringent, the effects are smaller. 

In the future, the observation of CP effects (or lack of them) in the — > v T oscillation can give new information on 
the light-heavy neutrino mixing. The bounds on such mixing, found in future neutrino experiments, can be stronger 
then the same bounds coming from charged lepton violating processes. 
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FIG. 1: Feynman diagrams for scattering of neutrinos in matter . All three diagrams contribute to neutrino-electron scattering 
m + e~ — > rife + e~ (/ = e). Only diagram (a) contributes to neutrino-nucleon scattering m + f — > rik + f (f = p,n). 
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FIG. 2: 

(a) The probability difference AP e <lf M defined in Eq.^ for L=250 km as a function of neutrino energy (oscillations in vacuum). 
The hatched region describes normal unitary oscillations with < 8 < 2ir. The shaded region corresponds to the non-unitary 
neutrino oscillations with additional CP breaking phases XiiX2 (Eq. ^) that change in the domain < \i 5i 27r. The e a 
parameters for both sets from Eqs. plp2| give e e e^ — 10~ 4 and the results are the same for both cases. Other oscillation 
parameters correspond to the best LMA fit values given in |2Ej (see the text for details). 

(b) Neutrino oscillations in matter with the same parameters as in Fig. 2a. In this case neutrinos pass only the first shell of 
the Earth's interior |24| with a constant density p = 2.6 g/cm 3 and Y e = 0.494. 

(c) Neutrino oscillations in vacuum with the same parameters as in Fig. 2a but for O13 = 0. In this case, in vacuum AP a ^n = 
for normal unitary neutrino oscillations. Non-zero effects are exclusively related to the non-decoupling heavy neutrinos. 

(d) AP^S^ for neutrino oscillations in matter for O13 — 0, L=732 km. Other parameters as in Fig. 2a. Because of matter 
effects, AP e % # even for unitary neutrino oscillations. 
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FIG. 3: APuS j ( a) and yl^_» T (b) for neutrino oscillations in vacuum, L—250 km. Other parameters as in Fig. 2a, e, 
according to Eq. pl[ 
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a) and A^^, T (b) for neutrino oscillations in vacuum, L=250 km. Other parameters as in Fig. 2a, e a 
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FIG. 5: APf?S T ( a) and A^!^ T (b) for neutrino oscillations in vacuum, L=732 km. Other parameters as in Fig. 2a, e a 
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FIG. 6: AP^St (a), A^H, T (b), and Aj^^ T (c) for neutrino oscillations in matter, L=732 km. Other parameters as in Fig. 2a, 
e a according to Eq. pTl 



